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jrt . Abtract: In this paper, we study a system of thermoelasticity with a degenerated 

second order operator in the Heat equation. We analyze the evolution of the energy 
density of a family of solutions. We consider two cases: when the set of points where 
^^D ' the ellipticity of the Heat operator fails is included in a hypersurface and when it 

^^ , is an open set. In the first case and under special assumptions, we prove that the 

evolution of the energy density is the one of a damped wave equation: propagation 
P I , along the rays of geometric optic and damping according to a microlocal process. In 

•^^ ' the second case, we show that the energy density propagates along rays which are 

• , distortions of the rays of geometric optic. 

-)— » 

C^ ' 1. Introduction 

We consider Vt an open subset of R and the following system of thermoelasticity 

d^u - Au + V • {i{x)6) ^ 0, (i, x) e R+ x fi, 
■ dte-V.{B{x)Ve)+-f{x).Vdtu^{), (t, x) e R+ X f7, 

^ ! (1-1) I W|t=o = "0 and dtu\t=Q = ui in fi 

^ ! [ W|an = 0, dtU\QQ, = and O^q^ = 0- 

i \ where u and 9 are scalar real- valued functions. The matrix- valued function x ^-^ 

^y ■ B{x) and the vector- valued function x i~7> 7(0;) are supposed to be defined on R** 

^~; \ (and thus on fi) and to depend smoothly on the variable x € R**. The matrix B{x) 

is assumed to be symmetric, non-negative: there exists C2 > such that 

(1.2) V (:e, e f^ X R'^, < B(a;)eC < C^2|Cr 

k>i ' Note that we may have det B{x) = 0. 

H 

Cu ' System (jl.ip has an energy 

E{u,e,t):= I \dtu{t,x)\^dx+ / \\Iu{t,x)\^ dx + / 9'^{t,x)dx 
Jn Jn Jn 

which decreases in time according to 

(1.3) E{u, 9, t) - E{u, 9,0)^-2 f f B{x)S/9{s, x).S/9{s, x)dsdx < 0. 

^0 Jn 

Equation (II. 3p gives a priori estimates for initial data uq G Hq(Q,), ui e L^(J7) and 
6*0 G i^(il) and yields by classical arguments the existence of a unique solution 
{u,9) e C"{Il+,H^{n)) nC\Ti+,L^{n)) x C°{R+,L^in)). We are interested in 
characterizing the way the energy decay; our aim is to describe the weak limits of 
the energy densities 

enit,x) = \dtUn{t,x)\^ + \W^Un{t,x)\^ + {9n{t,x)f, UEN, 
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2 A. ATALLAH AND C. FERMANIAN 

associated with families of solutions (Mn,,^n)neN of (jl.ip corresponding to families 
of initial data (uo,n)neN in one hand and (ui,n)neN, (^o,n)neN in the other hand, 
uniformly bounded in H^{^) and L^(r2) respectively. Without loss of generality, we 
suppose that (u")„gN and (0q )„g]sf goes to zero weakly in L^(0) and that (uo)„gN 
goes to zero weakly in H^{il). 

Main results on the subject are devoted to the situation where the matrix B 
is positive. It is known since the work of Dafermos [6 that the energy decays 
and the description of this decay has been the subject of several contributions. In 
particular, in [lS\ . G. Lebeau and E. Zuazua ^18^ have proved that the decay of 
solutions to ()l.ip is not uniform. They crucially use a result of Henry, Perissinoto 
and Lopes [15] which show that the semigroup associated to ()l.ip is equal up to 
a compact operator to the semi-group of a system consisting of a damped wave 
equation coupled with the equation on the temperature 9: 

, . f dfu - Am + Tu = 0, (i, x) G R+ x f7 

^ ' \ dt9-V ■{B{x)Ve)+-i{x).Vdtu^Q. 

The operator F is a damping operator given by 

r = G*Q^^G, G = 7(x) • V and Q = V • (S(a;)V-) . 

The behaviour of the energy density |M„(i,x)p associated with families of solutions 
of this damped wave equation (for initial data (wo,„)ngN and (wi^„)„gN uniformly 
bounded in H^{rt) and L^(51) respectively) can be studied in the same manner than 
in the paper of G. Lebeau and N. Burq (see [TB], [T7] and [J). In these papers, the 
damping is somehow different than the one of (|1.4p . however their method applies. 
One obtains that the energy propagates along the rays of geometric optic associated 
with the wave operator d^ — A with a damping depending simultaneously on the 
position and the speed of the trajectory. These trajectories are trajectories of the 
phase space {x(t),uj(t)) and the damping is microlocal: it is given by the principal 

symbol of the operator F, the function {x,uj) h^. — ^/^I'^l^ ■ 

The method of [1] is based on the use of microlocal defect measures in the 
spirit of the article [13] by P. Gerard and E. Leichtnam who are at the origin of 
that method. Then, similar works have been achieved for the Lame system in [4], 
for the equations of magnetoelasticity in [7] (see also [8]) and for the equation of 
viscoelastic waves by the authors [2]. 

Of course, the strategy that we have just described, fails if the kernel of B does 
not reduce to zero and we are interested in this situation: we assume that the set 

A = {ix,uj) eflx S'^"\ B{x)uj = 0} 

is not empty. Then, the operator V • (i?(x)V-) is no longer elliptic and something 
else has to be done. The method we use to treat the coupling between the tem- 
perature 9 and the amplitude u is mainly inspired by the analysis of semiclassical 
systems performed in |14) . Of course, we recover the result sketched in the previous 
paragraph when A = and we are also able to extend it to situations where A 7^ 
provided a weak degeneracy assumption stated below (see Assumption 12. 2p . This 
assumption consists first in a geometric assumption: the projection of A on J7 is 
included in a hypersurface T, and for (x, w) G A, the vector a; is transverse to S at 



ANALYSIS OF THE ENERGY DECAY OF A DEGEN. THERMOELASTICITY SYSTEM 



the point x. Then, Assumption 12.21 contains a compatibihty relation between the 
vector 7(0;) and the matrix B{x): 7(2;) G Ran_B(a;). With these assumptions, we 
are able to prove that the energy density is still damped along the rays of geometric 
optic even though they pass through E. At the opposite, if B{x) = in an open 
subset rj of rj, then the damping disappears and we have transport of the energy 
along rays which are distorsion of the rays observed before. Precise statements of 
our results are given in Section 2 and the organization of the paper is discussed at 
the end of this section. 

Notations. We will say that a sequence (un)„gN is u.b. in the functional space F if 
the sequence (un)nGN is a uniformly bounded family of F. We denote by |Ar|cd+2 
the hermitian norm of A" := (Ai, • • • , Xd+2) G C^+^i 

|A|(-;d+2 = A^ H h A^_|_2. 

Similarly, we will use the notation (AjF) for the hermitian scalar product of C^^^: 
VA = (Ai, • • • , Ad+2) € C^^', Vy = (n, • • • ,Yd+2) e c^'+'„ 

(A|y)c.+2 = XiTi + - ■ ■ + Xd+2Yd+2- 

2. Main results 

In this section, we present our results which crucially rely on the use of microlocal 
defect measures that we define in the first subsection: the evolution of the energy 
density is a corollary of the analysis of the behavior of microlocal defect measures 
associated with the sequence (un,0n)nGN- The second subsection is devoted to the 
analysis of properties of the thermoelasticity operator that are important for our 
purpose. Then, in the third subsection, we mainly consider the situation where 
the determinant of B vanishes on points of fJ which are simultaneously included 
in a hypersurface and in a compact subset of 51 (so that B is non negative in a 
neighborhood of 9f2). Finally, in the fourth subsection, we discuss what happens if 
B vanishes in an open subset of 57. 

2.1. Microlocal defect measures. Microlocal defect measures allow to treat qua- 
dratic quantities like the energy density by taking into account microlocal effects. 
They describe up to a subsequence the limit of quantities of the form {a{x, D)fm fn) 
where a{x,D) is a pseudodifferential operator and (/n)nGN a u.b. family of L^{fl) 
(or, more generally, of H''{n)). Recall that the pseudodifferential operator a(x, D) 
is characterized by its symbol a(a;, S,) which is a smooth function taken, for example, 
in the space A"^ of symbols of order m: this set contains the functions a = a{x,S^) 
of C°°{^ X R'') such that a is compactly supported in J7 as a function of x and 
satisfies 

Va,/?eN^ 3C„,^>0, V(x,0ef7xR^ |9^9f a(x,o| < C„,^(l + 1^1)""""- 

Then, the operator a{x, D), defined in the Weyl quantization by 

V/ G L^iR"), a{x,D)f{x) = {27t)-'' J a (^,^) f{y)e'i<^~y'>dyd^ 

maps 7J*(J7) into i7''^'"(51) (see [T]). We will also assume that for a e A"^, there 
exists a function aoo(a;,^) homogeneous of degree such that for 

(2.1) yiueS'^-\ Mm R-"'a{x,Ruj) = aoo{x,uj). 
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The symbols of Af^ are called interior symbols because they are compactly sup- 
ported inside fl. Such symbols are of no help for studying the behaviour of (/n)neN 
close to dfl: one then use tangential symbols which are defined in the Appendix. 
It is easy to convince oneself that the limits of quadratic quantities (a(x, D)fn, fn) 
(for (/n)rieN u-b. in L'^{fl) and a G A'^) only depends on the function Aoo. Then, 
following [12^ and ,20 , it is possible to prove that these limits are characterized 
by a positive Radon matrix- valued measure /x on f2 x S''^^ such that, up to the 
extraction of a subsequence, 

(a(x,D)/„, /„) — > (oocM)- 
Such a measure /i is called a microlocal defect measure of the family (/n)neN- 

Let us come back to the sequences (9tu„)„gN, (Va;'u„)„gN and (6'„)„gN that we 
need to study simultaneously. In that purpose, we set 

(dtUn 
On 

and the energy density e„ writes 

Gn[t,X) = \LJn(t,X)\f-,d+2- 

We shall consider quadratic quantities {a(x, D)Un{t), f7„(t)) where the symbol 

a{x,0 = ia^,j{x,D))^^^ 

is a {d + 2) X (d + 2) matrix of interior symbols of A™. Then, a microlocal defect 
measure of (C/n(i))nGN will be matrix- valued and will describe up to the extraction 
of a subsequence, the limits of the quantities {a{x,D)Un{t) , Un{t)). Of course, 
the i-dependence of these measures is an issue by itself. Therefore, since [/„ S 
L^Q^(R, L^(r2)), we test {a{x,D)Un{t),Un{t)) against smooth compactly supported 
functions of the variable t and consider the limits of 

I{a,e) := f e{t)ia{x,D)Unit),Unit))dt. 

A microlocal defect measure M of (t/n(0)neN is a positive Radon measure on the 
set R"*" X fl X S'^^^ such that, up to a subsequence, we have: Vd S C|^(R) 

Iia,e) — > / e{t)tr{a^{x,u})M{dt,dx,duj)). 

In particular, ¥6* eC5^(i?), VxeC;^(17), 



6{t)x{x)en{t,x)dxdt — > / e{t)x{x)tv M{dt,dx,duj). 

The matrix- valued measure -M(i, x, lj) is positive in the sense that its diagonal 
componants rriii are positive Radon measure and its off-diagonal componants rriij 
are absolutely continuous with respect to rrii^i and to^j. Using the special form of 
the componants of the vector C/„ (i) , one can write 

(mi{t,x,uj) mo^i{t,x,uj)uj mo^2{t,x,uj) 

iTio.i{t,x,ujYuj mo{t,x,uj)u! (^ oj mi^2{t,x,uj)u! 
rn^{t,x,uj) rrrr^(t,x,ujYuj vo(t,x,u;) 

where 
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• mi(t, X, uj), voit, X, Lu) are the microlocal defect measures of {dtUn)ne'N and 
of (0n)neN respectively, 

• mo{t,x,uj)u!iU!j is the joint measure of {dxiUn)ne'N and {dxjUn)nGTSi, 

• 1710^1(1, x,u!)ujj is the joint measure of (9tw„)„gN and {dxjU„)nG'N, 

• "T-o,2(i,a;,a;) is the joint measure of (6'„)„£n and (9tu„)„£N, 

• 'm'i.2{t,x,uj)ujj is the joint measure of (0n)neN and {dxjUn)neN- 

By "joint measure" of two u.b. famihes of i^(f^), (/„)tign and (5„)„gN, we mean 
a measure which describes the hmit up to extraction of a subsequence of quantities 
{a{x, D)fm gn)ne'N for symbols a of ordre 0. In the following, we assume that 
(t/„(0))„gN has only one microlocal defect measure M(0,x,a;). 

2.2. Analysis of the thermoelasticity operator. Let us now come back to our 
system (|l.ip that we rewrite as 

idtUn^Pix,D)Un 

where Un is defined in (12. 2p and 






«*V -iV • (7-) 


iV 





i-f-V 


iV-{B{x)V-) 



P{x,D) = 



We first study the eigenspaces of the matrix P(x, ^) which is not self adjoint if 
-B 7^ 0. However, we have the following proposition: 

Proposition 2.1. There exists Rq > such that for |^| > Rq, the matrix P(x,^) 
diagonalize with smooth eigenvalues and sm,ooth eigenprojectors. Besides, the kernel 
of P{x, ^) is of dimension greater or equal to d — 1. 

Denote by A_ , Aq and A+ the three eigenvalues of P{x, £,) which are not identically 0, 
and by 11", Xl" and n"*" the smooth associated eigenprojectors. Then, for all R > Rq 
and for x e C°°(R'') such that x(0 = for \l\ < 1/2, xiO = 1 V ICI > 1 and 
< X ^ 1; we have in 2?' (51), 



U>'ix,D)x{^]u,,{t,x 



+ o(l). 



(2.3) \U„it,x)\l,^.= Y. 

ke{o,~,+} 

This proposition is a consequence of Propositions 13.11 and 13.21 below where we 

study the asymptotics of P{x, f ) for large ^ (which differ whether I a;, -[It J G A or 

not). Each eigenprojector Il'^{x,£,) {k = 0,—,+) of P{x,(,) characterizes a mode 
and for each mode, we will analyze the microlocal defect measure of the componant 

(n'=(a;,i?){/„)„eN. 

2.3. Propagation and damping for weakly degenerated {B,j). Let us pre- 
cise first our assumptions. 

Assumption 2.2. We say that the pair (-6,7) is weakly degenerated if B and 7 
satisfy the following conditions. 

(1) There exists a hypersurface S o/R such that {deti?(x) = 0} C S and for 
all {x,U!) £ A, then uj is transverse to E in x. 

(2) There exists 7 G C°°(fl, R'*) such that 7(2;) = B{x)^{x). 



A. ATALLAH AND C. FERMANIAN 



Example 2.3. Suppose d = 2, Q = R, ^{x) = €2, B{x) = ( |-, 1 I • Suppose 
that b{y) — if and only if y — 0. Then, 

A-{(0,2/),(±l,0)),2/eR} 
and one can check that (B, 7) is weakly degenerated. 
The sequences 



"n " ~^^ I """:! ^ TTTT ' '^ x^r, 



1 / D 

have only one microlocal defect measure /i^ given by 

/iQ (a;,iuj) = mi(0, a;,aj) + mo(0,a;, w) ± Re (mo,i(0, x,uj)) . 

Theorem 2.4. Suppose that (-6,7) is weakly degenerated, that Supp(/XQ ) C A'^ 
and that there exists tq G R*+ such that for {x,uj) € Supp(/iQ ) and for t G [0, tq], 
{x±tuj,uj) e n. Then, for x € C(f([0,ro]) anrf G Cg°({detB(a;) 7^ 0}), 

hm / x{'t)4'{x)Gn{t,x)dtdx — / x{t)'t>{x)dlJ'f {x,uj)dt 



x{t)4>{x)dfi^ {x,uj)dt 

R+xOxS''-i 



w/iere /or a/Z a G C^f (f2 x S'*"!), 



(2.4) {a,^.t) = f aix±t.,.)E.^ - f iji-±^!^,. 



dfiQ{x,uj). 



Remark 2.5. 1) Note that even though the support of ip does not intersect the set 
{det_B(a;) = 0}, the trajectories x ±tuj which reach the support of (j) for t £ [0, tq] 
may pass through it. 

2) Besides, because of (1) in Assumvtion \2.Si a trajectory x + slo crosses T, at a 
finite number of times < ti < ■ ■ ■ < tpf < tq. Assuming B[x + tjUj)uj = 0, by (2) 
of Assumvtion [KR for j G {1, • • • , N}, there exists Cj G R such that 

(7(0; + SUj) ■ LU) . , 

-—, 7 — ^cAs — tA as s^tj. 

UJ ■ B{x + suj)uj ^^ '' ' 

This implies that the integral in {2.4^ is well defined. 



This result is proved in Section |42] The measures ^^ contain the part of the en- 
ergy corresponding to the projection of Un on the ±-mode. There is no contribution 
of the 0-mode (which corresponds to the temperature) because of the smoothing 
effect of the Heat equation. Note finally that the damping in (|2.4I) can be if for all 
times t G [0, tq], we have "f{x =F iw) • w = 0. We refer to section B^ for a discussion 
of what happens when (2) fails in Assumption 12.21 all the energy may be damped 
in finite time (see Remark l4.3p . 

Finally, in the Appendix, we shortly discuss what happens close to the boundary 
under the following assumptions: 

Assumption 2.6. (1) The rays of geometric optics have no contact of infinite 

order with the tangent to dQ. 
(2) There exists a compact K such that {deiB{x) = 0} C iiT C il. 
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(3) For all x e dVt, 7(2;) G T:,{dVt) 



If Assumption 12.61 holds, then one can use the methods developed in [13] and 
the papers [3] , [4] or [16] for the analysis at the boundary of microlocal measures 
of a family of solutions of a wave equation. One can prove the propagation of the 
energy along the generalized bicharacteristics as defined in [19] . We shortly explain 
in Appendix |X] how to reduce to the analysis of [1] , without writing in details the 
arguments which are exactly the same ones. 

2.4. Distorted propagation in an open set included in {B{x) = 0}. We 
suppose now that B{x) = in fJ an open subset of fi. Then, the symbol P{x, ^) is 
self-adjoint on 17 x R'' and the method of [TJ] can be adapted with straightforward 
modifications. We set 

and we consider fiQ and Dq the microlocal defect measures of the sequences 

^± = -^ (±u„4 - iW{x, D)\D\unfi + N{x, D)e^fl) 

9n - iN{x, D)\D\un,o + W{x, 0)6^,0 
where W{x, D) and N{x, D) are pseudodifFerential operators of order of symbols 

W{x,0 ^^^^x{^lRo) and N{x,0 = ^7^x(^/i?o); 

the function x is smooth and satisfies xiO = foi' ICI ^ 2 and x{C) = 1 ^^i \^\ > 4 
and Rq is chosen as in Proposition 12. II Note that the measures /Iq and i/q do not 
depend on the cut-off function x, besides, they satisfy 





'^o = „^^ , .^2 I M{x,uj) I (7(0;) • Uj) UJ 



c{x,L0) 
1 

2c{x,uj) 




C<i + 2 



^0 = ^ / „ I M{x,uj) I -UJ 



^±C{X,LU)^ 

—UJ 
j{x) -UJ 



Cd + 2 



'y{x) ■ UJ 

Before stating the result, let us introduce a notation. Define 

H,{x, = y^cix, e) • V, - V,c(a;, ^ • Vj, 

the Hamiltonian vector field associated with c and H^ the vector field induced by 
He on S*n. 

Theorem 2.7. Suppose that Supp(i/o) C 12, Supp(/Ig ) C 12 and that there exists 
To G R*+ such that for all x E fl and oj G S''^^, the projection on fl of the integral 
curve of H^ issued from {x,uj) stays in Q on the time interval [0,tq]. Then, for 
XGC(f(R) and(t)eC^in), 

lim / x{t)4'{x)cn{t,x)dtdx = I x(t)0(a;)(ipi(^(a;. ojjdi 

"^+°°jR+xn JR+xOxS'i-i 

X{t)(l){x)d'il:[ {x,uj)dt 

R+xOxSti-i 

x{t)(l){x)Ai'o{x,uj)dt 

R+xOxS<*-i 
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where ptj satisfy the transport equations dtjit T H^{x,uj)Ji^ = with initial 

1^0 



data ■■ 



Remark 2.8. Let us call distorted bicharacteristic curves the trajectories of S*fl 
associated with H^ : the transport equation for Ji implies that the energy propa- 
gates along these trajectories. 

Note that inside il, one cannot separate the contribution to the energy density 
of {un) and of (0„). The measures ptj depend of the value at time i = of both 
quantities (u„) and (0„). 

These results call for further works: it would be interesting to know what happens 
at the boundary of J7 and how transitions occur between the two regimes. It would 
be also interesting to know whether this result in {B{x) = 0} extend to the set A 
provided A is invariant by the distorted rays. The following example show such a 
situation. 

Example 2.9. Suppose d — 2, Vl — ^, 7(2;) — ei, B[x) — \ \^ 1 I • Suppose 
that b{y) — if and only if y = 0. Then, 

A^{{x,Lu), 3yeR, x^{y,0), l. = (±1,0)} 
is invariant by the distorted bicharacteristic curves issued from from points of A: 

{xf,^t)^{{y±sV2,0)A±hO)), seR. 

Note that, in that case, the distorted trajectories issued from points of A coincide 
with the usual ones; however, they are described with different speed. Note also 
that (2) of Assumvtions \2.2\ is not satisfied here. 

2.5. Organization of the paper. The main part of the article consists in the anal- 
ysis of the m.d.m.s associated with the families (U''{x,D)Un) (for k G {0,+,—}) 
where the functions n*^ are defined in Proposition 12.11 We begin in Section [3] by 
studying the symbol P{x,^), which allows to prove Proposition [2TTJ Then, in Sec- 
tion 01 we prove Theorems 12 .41 and 12 . 71 they rely on the analysis of the propagation 
of the m.d.m.s associated with the sequences (n^(x,D)Un) which is the object of 
Section [51 Finally, the Appendix [K\ is devoted to a discussion of the reflexion of the 
measures on the boundary. 

3. Analysis of the symbol of P{x, D) 

In this section we analyze the properties of the matrix P{x,(^). The main interest 
of Weyl quantization is that the symbol of a self-adjoint operator is real- valued. We 
denote by cr{A) the symbol of an operator A and we have in particular 





a(V) - i^, 


(3.1) 


(7(7(2;) -V) = il{x)-i--'^ -lix) 


(3.2) 


a{V-{B{x)V)) = -^B{x)i-^ + bo{x), 


(3.3) 


bo{x) = -\ i. dl,..B,,{x). 



l<j,k<a 
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Observe that, if rf = 1, the function bo has a sign on A. Indeed, if rf = 1, the 
points {x,uj) G A correspond to values x which are minima of B(x) and in this 
case 60 < 0. However, in higher dimension, the function 60(2;) can be positive or 
negative indifferently as shows the following example: choose d = 2 and B{x) such 
that A = {((0,y),(0,±l), y e R} with 

Bix)^[^^ ,^,.yOi\x\^) 

close to (0, 0). Then we have 

bo{x) = -^{12x1 + 6x2) +0{\x\^). 

Therefore, &o(0, y) = —^y + 0{\y\'^) and 60 changes of sign on A. Note also that if 

B{x) = in an open subset fl of f2, then bo = in f2. 

The eigenvalues of the matrix P{x, ^) satisfies the following proposition: 

Proposition 3.1. The matrix P{x,^) has a kernel of dimension greater or equal 
to d—1. Moreover, P{x,^) have three smooth eigenvalues A_, Aq and A+ with : 

1) // (x, ||t) ^ A, then dim [Kcr P(x, ^)] — d — 1 and 

(3.4) Ao(x,C) = -^^■B{x)C + ^bo{x)+^^^^^^^ + 0{\^\-'), 

with 

(3.5) /3(x,e) = |e|+0(l), ^(^^^)--2 j.B(^)^ +Om-'). 

2) // (a;, ill) e A, then, if bo{x) = 0, dim [KerP(x,C)] = d and if bo{x) ^ 0, 
dim [Ker P{x, ^)] — d ^ 1. Moreover, 

(3.6) Ao(a;,0 = ^bo{x)-^ + Om-'), 

X±ix,C) = ±l3{x,^)+ia{x,C) 
with 

and c{x,0 = ^/(^WW^W■ 

The modes ± (corresponding to the eigenvalues A±) give the wave feature of 
the equation. The speed of propagation is characterized by the function /? and the 
function a corresponds to the damping. Note that in the second case, the speed 
of propagation are distorted by comparison with the initial wave operator d^ — A. 
Outside A, the eigenvalue Aq encounters of the Heat aspect. 



(3.7) /3(x,e)-c(a;, 0+0(1), a(x,g) ^ -&o(x) ^ ^^J \^ +0(1^1'^), 



Proof We write P{x,^) — iQ{x,^) and for simplicity we work with Q{x,^). For 
p,q E N*, we denote by Op.^ the p x q matrix with all coefficients equal to 0. We 

have 

i^ -k{x,0^ 



Q(x,e)= z*e Od,d 



d.l 



,k{x,0 Oi,d bix,0 
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where, in view of p.ip and p.2p . 6 = 62 + 5o with 62 — —B{x)£_ ■ ^ real-valued and 
k = iki + fco with 

ki^-l{x)-£,, fco = -V-7(a:). 

The vector {x, Y,y) e R x R'' x R is an eigenvector of Q for the eigenvalue 1/ if and 
only if 

(1) lyx - i(,-Y + ky = 0, 

(2) ix^ - vY =0, 

(3) kx + {b-v)y = 0. 

Let us suppose first that v — Q, then for ^ 7^ 0, equation (2) gives a; = 0. For 
& ^ 0, equation (3) then implies y — Q and (1) gives ^ • F = 0. Therefore, for large ^ 
the kernel of Q{x, f ) is of dimension at least d—\ and contains the smooth vectors 

^.■(e) = (0,e,(O,0) 

where {s.j{S,))i<j<d is a smooth basis of the hyperplane orthogonal to ^ (such a 
basis exists for ^ 7^ 0). 

Let us now suppose v ^ Q. Equation (2) implies that Y is colinear to ^ and 
X = —ivr where Y = r^. Equations (1) and (3) become 

ky - i{i^^ + |$|2)r = 0, 

{b — i')y — iukr = 0. 

Therefore, the non zero eigenvalues are the roots of the real-valued polynomial 

f{X) = ~X^ + X% - X{\e -f Ifcp) + 6|^|2. 

It is easy to see that for large ^, / has one real- valued root v^ (with vq ^ 
for 60 7^ 0) and two conjugated complex- valued roots v+ and v^. These three 
roots of the polynomial / are simple and thus smooth; they give three smooth 
eigenvalues of Q. Consider three associated eigenvectors Vb, V+ and V-, they are 
independent from the vectors Vj{S,) defined above. Therefore, we are left with a 
basis of eigenvectors: the matrix Q{x^S,) diagonalize. 

Let us now study more precisely the asymptotics of the eigenvalues. 
1) Suppose (x,^) ^ A, ^ large and denote by Xi < X2 < the two negative roots 
of f'{X). We have Xi = §62 -f 0(|C|) and X2 = 0{\i\). Then, using /'(X2) = 0, 
we obtain /(X2) = 6(X| -I- |^p) < and deduce that / has only one real- valued 
root vq with vq < Xi. We set vq = cf)b with </) > | + o(l) and 

(3.8) = /(z^o) = b^4>\l - 0) - 06(|fc|' + len + 6|er 

Necessarily, 4> = l-fr withr = -kl/bl + O {\i\-^) , whence i/q = b-^ + 0{\i\-^) 

and ([Ol . 

Let v± be the two other (non real-valued) roots, we set v± — a^ i/S. We observe 

i^+ + i^_ = 2a = 6 — i/Q, 



whence a = §- + 0{\^\-^) and /S^ = |^|2 + 0(|^|). This implies ^ 
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2) Suppose now (x, f ) E A, then b{x, £,) — bo{x). The polynomial function ,f'{X) has 
no real- valued root and f{X) has only one real- valued root vq. Since /(0)/(feo) < 0, 
the function 1/0(2:,^) = 0(1) as ^ grows and 

M\^f + \kf) = bo\e + 0{l), 

whence (|3.6p . Besides, denoting as before by i/± = a^i/3 the two other roots, p.9p 
gives 

" = I fcfTliP + ^^'^'"'^ ^""^ /?' = A:? -H I^P + 0(1), 
whence (|3.7I) . <^ 

Let us now describe the eigenspaces of P{x,^) for large ^. 

Proposition 3.2. For ^ large enough, there exists a smooth basis of C'^'^^'^ , 

{Vo,V+,V-,Vj), l<j<d-l such that 

P(x,0Vfe = Afel4 for e{0,+,-}, 

P{x,OVj=0 for l<j<d-l, 

with the following expansions 

1) // (x, ^) ^ A, 

(3.10) IS = (O,ej(0,0), 6,(0 -e-O, l<j<d-l, 

(3.11) y± = _!= (^^1,^,0)+ o(ier^), 

(3.12) Fo = (o,o,i)+o(ier'). 

2) // (a;, ^) e A, 

(3.13) IS- = (0,6,(6,0), e,(e)-e = 0, l<j<d-l, 

V2 V c(a;,0 c(a;,6y 



(3.14) 



(3.15) Vo = (o,^^,^)+0(|^ri). 

Note that there exits smooth eigenvectors but their asmptotics are discontinuous; 
similarly, even though they are not orthogonal projectors, their asymptotics give a 
decomposition of C^^^^^ on orthogonal subspaces. Moreover, when bo = 0, there is 
an eigenvalue crossing between the eigenspace for Aq which merges into the kernel 
of P. However, there still exists a smooth basis of eigenvectors and we will take 
advantage of this fact in the following sections. Besides, we have the following 
remark. 

Remark 3.3. Assuming (2) of Assumvtion W^ we have for {x,uj) £ A, 

j{x) ■ U! = ^{x) ■ {B{x)u!) — 0. 
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The 



efore, if [x, j|| j G A, c{x,^) = \^\ 



Vj = (0,6,(0,0), e,(O-e-0, l<j<d-l, 



y± = 4f-i'±JT'0)+o(i^r'), 



Fo = (o,o,i)+o(icr^) 

Let us now prove Proposition [3?2l 



Proof : In the proof of Proposition 13. 1[ we have already buih the vectors Vj . 
Besides, we have seen that the eigenvectors of Q{x^^) associated with vq, vj^ and 
i/_ are of the form [x, rS,, y) with x = —ivr and 

(6 — i>)y — ivkr = 0, 

ky - i(i/2 + |^|2)r = 0. 

Let us consider first the ±-modes. We have i'± = 0(|CI)- Therefore, 6 — i/± ^ 
for large ^ independently of the fact that (a;,C/|C|) G A or not. The vectors 

V±{x,£,) = r{iy±{iy± - b),i{iy± - b)£„kiy±) , 

1/9 

with f = (|j/± — fep(|j^±P + I^P) + |fcp|i^±P) are smooth non-zero eigenvec- 

tors associated with h'± . In view of the asymptotics of Proposition 13.11 we obtain 
asymptotics for V±: 

(1) In A^ i/± = T«|^| + 0(1), k = -ij{x) ■ e + 0(1) and b = ~B{x)C ■ C + 0(1). 
Therefore 

f^l=i\c\Bix)(-o-'ii + om-'))^ 

V±ix, = r [{T^^\{B{x)^ ■ 0,^iBix)i ■ C, T|ei7(a;) ' + 0(1^1')] , 

whence p. 111) . 

(2) InA, z/± =T'ic{x,^)+0{l), 6 = 0(1) and we still have fc = -ij{x)-^+0{l). 
Therefore 

f=l/iV2c{x,0'){l + 0{l)), 

V± {x, = r [{-c{x, C)\±c{x, C, Tc{x, Oli^) ' + 0(|C|)] , 

whence ([314]) . 

Let us consider now the 0-mode. We have t^g = 0(1) in A and t^o + 0(1) G R 
in A^. Therefore, i^g + |^p 7^ for large £_ and the vector 

yo(a:,0 = ^0 {-iyok,-ik^,iy^ + |CP) 
with ro = (l^oHfep + \e\k\^ + h' + 1^1')')"'^' 

is a smooth eigenvector associated with the eigenvalue i/q. We are now left with a 
smooth basis of eigenvectors. Let us now study the asymptotics of this vector. 

(1) In A^ 1^0 = -B{x)^ ■ C + 0(1), therefore 

ro = iB{x)^ ■ 0-' (1 + Oder ^)) and ^0(2;, = (0, 0, 1) + Od^p^). 
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(2) In A, 1/0 = 0(1), whence 

ro^i/{cix,o\mi+om-'), 

Vo{x,0=ro[iO,ili^)-Oi,\i\')+Om)], 
which gives p.lSp .O 



Before concluding this section, we point out that these asymptotics imply Propo- 
sition [2TlJ Indeed, we have obtained the existence of Rq such that for |^| > Rq, 
P{x, ^) diagonalize with smooth eigenprojectors given by the Gram matrix of the 
basis of eigenvectors of Proposition 13.21 Asymptotically, these projectors are or- 
thogonal and we have 

\i\x,0 = Vk{x,0®Vd^)+o{l). 

Since Un{t, x) goes weakly to in i^(r2), we have 

2 



\Un\t, x)\(2d+2 



XI -^)f/„(i,a:) 



o(l) in V'{^). 



Besides for R > i?o, we observe that Ii^{x,D)x{D/R)Un{tTx) is well-defined and 
we have in I?'(r2), 



Ii\x,D)x\^]Un{t,x) 



Cd + 2 



Vk{x,D) 



X ( -^ ) Un{t,x] 



Cd+2 



oil) 



We can now use the asymptotics of Proposition |33l which, combined with the weak 
convergence to of t/„(i, x) gives 

Y, {vu{x,D)\x(^]ur.{t,x)] 
fee{o,-.+} ^ ^^^ ^c^+^ 

= {dtu,,{t,x)f + \Vun{t,x)\' + {e^{t,x)f + o{l) in v'{n), 

whence Proposition 1 2. II 

4. Proof of the main results 

The proofs of Theorem 12.41 and Theorem 12.71 are inspired by the method de- 
veloped in |14| for analyzing semi-classical measures associated with solutions of a 
system of p.d.e's. The proof of Theorem 12.71 is a direct adaptation of the results 
of [15 in the microfocal defect measures setting while the proof of Theorem 12.41 
requires non-trivial adaptations due to the fact that P(x, ^) is not self-adjoint and 
that one of its eigenvalue is a symbol of order 2. Therefore, we focus on the proof 
of Theorem 12.41 and we let to the reader the simple adaptation of these arguments 
to prove Theorem 12.71 Theorem 12.41 relies on Propositions 14.11 and 14.21 stated in 
Subsection 14.11 then the proof of Theorem 12.41 is performed in Subsection 



4.1. Preliminaries. We state technical results (Propositions 14. l l and 14.21 ) that we 
will use in the next subsection. The proof of Proposition 14.11 is done at the end of 
this section while the one of Proposition 14.21 is postponed to Section [5] 



The first result describes the evolution of the temperature 0„ : 

Proposition 4.1. Let v be a microlocal defect measure of the sequence (0n)neN of 
Lf^^(R,L'^{'R'^)), then, as a measure on R+ x A'^, j/ = 0. 
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Our second result concerns the contribution to the energy density of the sequences 
{Il^{x, D)x{D/R)un)^^j^ for R > Rq and x a-s in Proposition l2.1l We set 

C/„%(i,x)=n±(a;,i?)x(^|)[/„. 

The sequences (t/^j /j)neN are uniformly bounded in L^{D,) and their microlocal 
defect measures are matrix- valued measures independent of i? > 0. Besides, by the 
definition of U^j^, these measures are of the form /x-|-(i,x,w)n (a;,a;) where the 
measures iJ,^{t,x,uj) can be understood as the traces of fi±{t,x,uj)Il^{x,uj). We 
prove the following result. 

Proposition 4.2. Assume condition (2) of Assumption [27B is satisfied in fli d fl 
and let T > 0. There exists a subsequence nt and a continuous map t i— )■ fJ-±{t) 
from [0,r] into the set of positive Radon measures on ft x S"*^^ such that for all 
t € [0, T] and for all scalar symbol a e A'^, we have 

(4.1) (a(a:,i^.)C/,t,flW|C^„liiW) ,— > f a^ix,uj)dfi±it,x,uj). 

Moreover, in V {{t > 0} x f^i x S**^^), we have 

(4.2) 9t/i± ± Hfs{^i±) - 2aooH± = i'± 

where v± is a measure supported on A absolutely continuous with respect to /i± and 
where for all a € J^, 

a^{x,U})d{Hi3{fj,±)) = - / {Hpa)aodii± = / ({a, /3})oo d^±. 

Proposition 14. 21 is proved in Section [5] below. 
Let us now prove Proposition 14. II 



Proof. It is enough to prove that if g is a symbol of order such that q{x,uj) G 



C°°(A^), then q{x, D)9n goes to in i^^^(R, L'^{Vt)). We observe that we only need 



to consider large values of ^. Indeed, if x £ C°°(R''), x(C) = for 151 ^ 1 and 
xiO — 1 for 1^1 > 2 with < x < Ij we have 

q{x, D)en = q{x, D)x (^ 9,, + o(l) 

in L^{n,) because 9n goes weakly to and the operator 1 — x ("§) i^ compact. We 
write 



q{x,Ox (I) = ^QB.ix,0 ■ VW)i 



where Qr is the vector- valued symbol of order —1: 

with x{'^) — W\~^x{'^)- Note that Qr is smooth since g = in a neighborhood 
of A and because x G Co°(R-'^)- Note also that Qr satisfies symbols estimates 
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uniformly with respect to R. Therefore, there exists a constant C > such that 
ioT all f e L^ifl), 

\\QRix,D)f\\mn) <C\\f Wmny 

We write 



qix, D)e,, - R-^Qb{x, D)^B(x)V^en + o(l) 
in L^(r2) where the o(l) is uniform in i? > 1 when n goes to +00. As a consequence, 
if / e L'^{n) and V e C^(R), we have 



V'(i)(<z(x,i^)^„(t)|/)di 
1 



7/^(0 /BMV0„(i) I Qnix^DYf] dt 



0(1) 



C 



<j^\\f\\mn)J |V(i)lll/B(^V0„(i)llL^(o)di + o(l) 
for some constant C > and where the o(l) is uniform in i? > 1 as n goes to +00. 



We observe that the energy equahty (jl.Sp gives that the family ( ^J B{x)V xOn{t) 

is u. b. in L^ ([0,T] x Vt,Il'^) for aU T > 0. Therefore, letting n and i? go to cx), 
we obtain the result. (^ 

4.2. Proof of Theorem 12.41 We suppose that (i?,7) is weakly degenerated, that 
is that the conditions (1), (2) and (3) of Assumption 12.21 are satisfied. Let X: 'P be 
as in Theorem 12.41 We observe that the energy density e„(i,a;) writes 



e„(t,a;) = |0„(t,a;)|2 



^ Y7 

- 7^ ■ VU„ 



StMn + i-TTj- • Vu^ 



Therefore, it remains to analyze the limit of each of these terms for (j) supported 
outside {det-B(a;) — 0}. 

First, we observe that by ProDOsition l4.1[ if (/> e C^ ({deti?(x) 7^ 0}), 

x{t)4'{x)dv{t,x,uj) =0. 



lim 

n— >+oo 



X{t)(l){x)\en{t,x)\'^dtdx 

iRJn JrJa 

Therefore, the weak limit of the energy density express only in terms of the se- 
quences dfUn ± -Mr ■ Vu„, of which microlocal defect measures are iJ,±{t,x,uj) by 
Remark 13.31 By Proposition 14.21 the measures /i± satisfy 

dtfJ-± ± wVx/i± = -Iac —- ti± + v± 

B{x)(jj ■ a; 

where we have used Q;oo(a;,w)lA = ^60 (■'^) efa Ji*^ = by Remark 13.31 
Let us now prove that the fact that lo is transverse to the hypersurface E, implies 
that /i±ls = 0. Let f{x) = be a local equation of E in a subset VI2 of Vi and let 
a be a symbol supported in 0,2 x R"* and such that a > and ^ • V/(a;) > for 
(a;,^) £ Suppa, ^ 7^ 0. We choose a function E Co(R) such that d'{0) = 1 and we 
use the test symbol 

bs{x,O^Saix,O0' ^^""^ 
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where S > 0. The transport equations for /i± imply 



a^{x,uj)uj-Vf{x)0' (^ ] dfi±{x,uj) = 0{S). 



By letting S go to 0, we obtain J,,.f^aoo{x,uj)uj ■ \'f{x)d^±{x,uj) = 0, whence 

/ij-ls = on the support of a. By this way (inspired from j9]), we finally obtain 

/i±1a = since A C S x R''. 

Besides, since the measure i'± is supported on A and absolutely continuous with 

respect to /i±, we deduce iy± = 0. 

Finally, we observe that the function F{x^lo) = — '^(J\^l^ extends continuously to 

the set fJxS'^-i with i^(a;,cj) =OonA (since 7(2:) -w = ^{x)-{B{x)ijj) = 0{\B{x)lo\) 

by (2) of Assumptions I2.2p . Therefore, we can write 

„ , T-, (7(2^) ■ ^Y 
^t^i± ± wV^Mi = — — — ^l±. 

B{x)(jj ■ (jj 

As a conclusion, we obtain Theorem 12.41 Indeed, take a E C^{il. x S''^^), < s < 
t < tq and set 

ft-S 

at{s, x,uj) = a{x ± (t — s)w, w) Exp 



pt — s 

/ F{x ± aijj,uj)da 
Jo 



then at(t, x, to) = a{x, w) 



whence formula p.4p . 

Let us conclude this section by a remark. 

Remark 4.3. // (2) of Assumption \2^ fails, then all the energy may be damped 
in finite time. 

Proof: Suppose /ip = and /i(j" = ao5{x — xq) ® (5(a; — wq) for {xq,ojq) ^ A and 
suppose that there exists tq such that 

Vt e [0, io[, (a;o + ii^o, i^o) ^ A and (xq + ioi^o, i^o) e A. 

Set yo = xq + tgujo and suppose that dB{yo)[uJa, ujq, luq] 7^ and 7(2/0) • ^^o 7^ then 

'■* (7(0:0 + sujo) ■ wo)^ 



cxp 



ds 



/o B{xo + sujo)ujQ ■ ujo 
and all the energy is damped between times and tg 



— > 



5. Propagation of Microlocal defect measures 



In this section, we prove Proposition l4.2l We consider the + mode; the — mode 
can be treated in the same way. We proceed in three steps: 

• First, we analyze the time derivative of (a(x,-D)[r^^(i) | U^ jj{t)] for scalar 
valued symbols a E A^ and prove that there exists a symbol T E Ai such 
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that, uniformly in i? as n goes to +cxi, 

(5-1) |(«(^>^)f/+fi|f/+«) - -{{'^^P}{^,D)U+JUI^) 

+2{{aa){x,D)Ul^\Ul^) 

T{x, D)x (^^) t/„ I X (^) U,}j + 0(1). 

• We calculate precisely the symbol T{x, ^) and show that T E A'^. Therefore, 

the quantity ( ^ [a{x, D)U^ ^ | ?/„ ^1 ) is uniformly bounded and by 

considering a dense subset of A\^^ Ascoli Theorem yields the existence of 
the continuous map i H> ^± [t) satisfying (|4.ip . 

• Finally, we prove that for all f/; G C(j"(R), 

(5.2) Jhn^ j i;{t)[T{x,D)x[^Ur.{t) \ x{^U.n{t)^dt 

tl;{t)aoo{x, uj)di'±{x, uj)dt 



where v± is a measure supported on A and absolutely continuous with 
respect to /i± 

At the end of these three steps, we have obtained Proposition 14.21 We now detail 
each of these steps. 

5.1. 1st. step: proof of (|5.ip . The family of functions U^ j^ is solution of the 
equation 



idtU+j,{t, x) = U+{x, D)x [^ P{x, D) 



(5.3) idtU+Jt, x) = n+ix, D)xi-] P{x, D)Un. 



Since Vx is compactly supported, we have 

X (^^) P{x, D) = P(x, D)x [^ + Kr{x, D), 

where Kj^ is a compact operator. Moreover we have 

n+{x,D)P{x,D) ^ X+{x, D)U+ {x, D) + Ri{x, D), 

where the symbol Ri will be precisely calculated in the next subsection. So ([5 
becomes 

tdtU+j,{t,x) = {X+{x,D)U+{x,D)+Ri{x,D))x(^]ur, 

+U+{x,D)KRix,D)Un 

= X+{x,D)U+j^ + Ri{x,D)x(^jUn + n+{x,D)KBix,D)Un 

(5.4) = X+ix,D)U+j, + F+j,{x,D), 

where 



(5.5) 



F+j,{x,D) = Ri{x,D)x(^]un + n+{x,D)KR{x,D)Un 
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For a real-valued symbol a e A'^, (|5.4p implies -ri^i^^^W^ rIU^ r) = h + h, 
with 

h - \{a{x,D)\+{x,D)Ulj,\Ul^)--M^,D)U+j,\\+{x,D)Ul^), 

h = \{a{x^D)Flj,{x,D)\Ulj,)--Mx,D)Ul^\Fl^). 

The term /i will give the transport by the vector field Hp and the damping by a. 
The term I2 is a rest term and its main contribution will be described by a symbol T . 

Let us study /i. 

h = ^{{a{x,D)X+{x,D)-{X+{x,D)ra{x,Dj)U+jU+j,). 

We recall that X+{x,D) — (3{x,D) +ia{x,D), so, since we use the Weyl quantifi- 
cation for the symbols we have (A+(a:, D))* = j3{x, D) — ia{x, D) and 

h = ^{[aix,D),l3{x,D)]U+jU+J+2{{aa)ix,D)U+jU+J 

+ {r^,{x,D)U+^\U+J 
= ~{{a,(3}ix,D)U+jU+J+2{{aa)ix,D)U+jU+J + (r_i(x,I?)f/+^|C/+«), 

where r_i(a;,^) will denote from now on a generic symbol in A^ . Using p.Sp . we 
then obtain if n, R tend to infinity 

(5.6) /i ^ - {a,(3}oo{x,uj)dfi+{x,uj) +2 / aoo(a;,i^)Q:oo(a;,a;)(i^+(a;,a;). 

Note that, for these terms, we do not need to integrate in time to get the conver- 
gence. 

Let us now study l2- We set Un,R — X (§) Un and bv l5.51 we obtain I2 — l2,i+h.2 
with 

l2,i[t) = -{a{x,D)I{+{x,D)KRUn{t)\Ii+{x,D)Un,R{t)) 
I 

--{a{x,D)I]+{x,D)U,,^R{t)\I]+{x,D)KRUAt)) 

I 

h.2(t) = -{aix,D)Riix,D)Un.R{t)\Tl+{x,D)Un.R{t)) 

--{aix,D)n+ix,D)Un,Rit)\Rii^.D)Ur,,Rit)) 

I 

Since (C/n)„gN goes weakly to as n goes to -l-oo and Kr is a compact operator, 
the sequence {KRUn)neT<! goes strongly to 0. Therefore 

(5.7) VteR+, l2.i{t)^0. 
Besides, h^it) - {T{x,D)U,,^r \ U^,r) with 

Tix.D) := - in+{x,D)*a{x,D)Ri{x,D)-Ri{x,D)*a{x,D)U+{x,D)) . 
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5.2. 2nd. step: the symbol of the rest term. In the following, rest terms in 
the symbol class A^ for fc e N will be denoted by r_fe. We prove the following 
Lemma. 

Lemma 5.1. Suppose f2i is like in Proposition \4-2\ Then, in fli, T G A^ and 

T{x,0 ^ Ti{x,0 +T2ix,0 + r-i{x,C) 
with 

(5.8) Ti - -^(n+{n+,/3} + {n+,/3}n+ + n+{n+,/3}n- 

+u-{n+,(3}u^ 

(5.9) T2 - ^{n+{n+,b}K-K{Il+,b}Il+) 

2,1 

Besides, TiHq and HqTi are symbols of A^ . 

Note that Tj G A*^ because 11+ K, Kn+ £ A'^^ by Remark 

Proof 

We first calculate R\ , recall that 

Ri{x, D) = n+(.T, D)P{x, D) - \+{x, D)Ii+{x, D). 
We write P{x, = Pii^^ + M^, CjK, with 

Pi{x,0=\ i 0,,, 0,,i andif=(0™+i ^Y'' 

Since b{x, ^) is of order 2, i?i is a priori of order 1 and in view of I\+P = A+n+, 
we have 

i?i = ^ ({n+, P} - {A+, n+}) + QK + r_i 

where the matrix- valued symbol Q is in A^ and is the sum of second derivatives 
of n+ multiplied by second derivatives of b. More precisely and after ordering the 
terms of higher degrees, we write 

Ri = ^{n+ ,b}K-'- ({n+,/3} + {n+,Pi}) + g/^ + r_i, 

where we have used that \^ = (5 + ia and {a, 11+} e -4^^^. However, since 7(2:) S 
Ran_B(a;), then 7(0;) • ^ = for (a;, -[It 1 G A and, by Remark 13.31 

(5.10) Kn+ ,n+ K e A~^ . 

Derivations of these relations imply {I\+ ,b}K e A^ and QK e A~^ . Finally, we 
find 

i?i = i{n+,6}i^- ^ ({n+,/3} + {n+,Pi}) +r_i e A°, 

whence T £ A^ with 

T = i(n+ai?i-i?5;an+) +r_i 

i ' 

= ^(n+{n+,fe}/^-x{n+,fe}n+) 

-^ (n+{n+,/3} + {n+,/3}n+ + n+{n+,Pi} - {Pi,n+}n+) +r_i. 
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where we have used (11+)* = 11+ + r_i and 

Rl ^ \k{iv-^M + \ ({n+,/3} - {Pi,n+}) + r_i. 

We now transform the expression of the principal symbol of T. We write 
Pi= P-ibK = A+n+ + A_n" + AqH" - ibK. 

By Proposition 13. 1[ 13.21 and Remark 13.31 we notice that 

• in A^ no = if + O(lCr^) and Xq = ib + 0(1), 

• in A, n° = K + O(ler^), Ao = 0(1) and b = (9(1). 

Therefore, 

Pi = A+n+ + A_n- + ro. 

In view of I^Mi we deduce T = T2 + T with 

f = -^(^(n+{n+,/3} + {n+,/3}n+) 

+ J2 (n+{n+,A,n^}-{A,tf,n+}n+)]+r_i 

= -^(^(n+{n+,/3} + {n+,/3}n+) 

+ ^ (^+{^+,/3^^}-{/3^^^+}^+)')+r_l. 

Observing that {n±,n+}, {n+, n±} e A~^ by Remark [331 and equation ([51]), we 
obtain 

T = -^(n+{n+,/3} + {n+,/3}n++ ^ (n+{n+,/3}tf -tf{/3,n+}n+) 

= Ti + r_i 

where we have used n+{/3,n+}n+ = (which comes from (11+)^ = 11+ whence 
{/3,n+} = n+{/?,n+} + {/3,n+}n+ and, multiplying by n+ on both sides, we 
obtain n+{;3,n+}n+ = 0). Notice that {/3,n+}n" G A'"^ and nO{/3,n+} e ^;"^ 
whence n°T, Tn° e A~\ 



5.3. 3rd. step: passing to the limit in the rest term. We use the following 
Lemma. 

Lemma 5.2. Consider A a smooth symbol of order supported in Qi and such 
that A = U+An^ or A = nMn+ where 0^=0" or H^ = H^ 1 < j < rf - 1 . Then 
z/VGCo°°(R), 

%(t)(^(:E,i5)C/„(i) , C/„(i))rft —> 0. 

n— >-+CX3 
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Since the symbol Ti is the sum of terms of the form n+{n+, /3}n^ or U'^{U+,l3}U+ 
with £ = — or i = j, I < j < d — 1, we can apply the Lemma and we obtain 

V'(t) (t{x, D)x (^I) Unit) I X (f ) Un{t)^ dt 

i:{t) (t2{x, D)x (^ Unit) I X (^) Un{t)^ dt + 0(1) 
as n goes to +00. In view of l5.9[ we have 

Ihnanp J ^/j{t) {t2{x,D)x [^ ^"W 1^ (f) ^"W) ^^ 



, x{t){T2)oo{x,uj)dv+{t,x,uj) 

where v.^ is the joint measure of KUn = ^n and of U^j^. The measure i/_|- is 
absolutely continous with respect to /i_|_ and h'. By Proposition 231 ^'+ is supported 
on A and we obtain (J5.2I) with i^-)_ = (72)00^'+ which is absolutely continuous with 
respect to ^+ and supported on A. We now focus on the proof of Lemma 15.21 



Proof 

Let n^ be one of the projectors 11^ or 11^, 1 < j < d — 1 and let us consider a term 
of the form H+AH^; the proof is similar for the other terms. We have 

[n+Arf , P] = (A^ - A+)n+^tf. 

Let us denote by C^ the symbol of A^^ 

C^ = {Xe-X+)-^n+AU'^. 

We have C"^ e A:[^ and, by Remark[331 we have C'^K, KC^ e A^"^ . We write 

j ^{t) {{n+ An'){x,D)UnAUn,R) i^.^^.^dt 

^ J m{{[C\P])ix,D)Un,B.\Un,I,)^.^^,^dt. 

Besides, 

[C^ P]{x, D) = [C\x, D), -idt + P(x, D)] - 1 ({C^ P} - {F, C^}) (x, D) 

+r^i{x,D). 
Using P = Pi + i6is: and {C^ Pi}, {Pi, C""} £ A^\ we obtain, 

{C^P} - {P,C'} = z{C^6}i^- zif{6,C^} + r_i. 
Therefore, we have 

[C',P]{x,D) = [C^(x,i^),-za,+P(x,I?)]-i({C^6}if-X{fe,C^})(x,i^) 

+r_i(a;,i:') 
= [C\x,D),-idt+Pix,D)]+r^i{x,D) 

where we have used {C'^,b}K = {C'^K,b} e A^^ and K{C\b} = {KC\b} e 
A-^ (as a consequence of KC\C^K e A:[^). In view of P*{x,D) = P(a;, Z?) - 
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2ib{x, D)K , we can write 

[C^ P]ix, D) = C\x, D) i-idt + P{x, D)) - i-idt + P*(.T, D)) C\x, D) 

-2ib{x, D)KC\x, D) + r_i(x, D). 
Finally, using [idt — P{x,D))Un,R = KjiUn with Kr compact, we obtain 

j X{t) {{^-^ An'){x,D)Un,R\Un,R) ^,^^,~^dt 

= -ijx'it) {C'{x,D)U^^R\U^^R)^,^^,^dt 

-2i J x{t) {h{x, D)KC'{x, D)Un,R\I<Un,R)^,^-^,.^ dt + o(l). 
Since C^ is of order lower or equal to — 1 , we have 

We use 

b{x, D)KUnM =(o,--- ,0,x (^] b{x, D)0n] 

with b[x,D)dn e Lfoc(R, if^^)) by Proposition I4l1 Since moreover bKC^ is of 
order in fJi, we have 

Yimsnx^ f ^{t){b{x,D)KC\x,D)U^,n\KU,,^R) dt -^ 0, 

whence 

limsup j ^{t) {{W+AU^){x,D)UnM\UnM)L2,^.)^i 



R~^+oo 

Appendix A. Analysis on the boundary 

In this appendix, our aim is to prove that, under Assumtion 12.61 the energy 
reflects on the boundary according to the laws of geometric optics. We shortly 
recall the arguments of [4] and explain how they apply to our setting. In all this 
section, we work in a neighborhood fii of a point of dVt. We first recall in the first 
subsection the definition of Melrose-Sjostrand compressed bundle (see [19] and the 
survey [3]) and of the generalized bicharacteristics. Then, in the second one, we 
will link [4]'s approach and ours (we also refer to [3] and [7]). The last subsection 
will be devoted to the proof of the main statement of this Appendix. 

A.l. Melrose- Sjostrandt compressed bundle and the generalized bichar- 
acteristics. We work in space-time variables and set L = Rj x Jl. We denote 
by (2:,C) the points of r*L: z = {t,x) and ( = (t, ^). Then, the Melrose-Sjostrandt 
compressed bundle to L is given by 

T;L ^ {T*L \ {0}) U {T*dL \ {0}) . 

Quotienting by the action of R^ through homotheties, one obtains the normal 
compressed bundle to J7 

s;l = t;l/r+. 

The projection 

vr : (T*R''+i)|-^->Tfc*L 
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induces a topology on Tj^L. On T*L, we denote by po the symbol of the wave 
operator and by So the projection on T^L of its characteristic set 

(A.l) Po{z,(:) = r'-\e, So = ^({(z,C), Po{z,C)^0}). 

Locally, near a point of dL, we use normal geodesic coordinates 

{y,v) = (yo,yi,--- ,yd;vo,vi,--- ,Vd) e 'R'^''+'^ 

in an open set V of r*R''+i such that LnV = {yd > 0} nV. Therefore, the wave 
operator d^—A writes —dy^ — R{y, Dyi) where the symbol R{y, ij') is a homogeneous 
polynomial of degree 2 in 77'; we denote by r{y, -q') its principal symbol. We can now 
distinguish between different sorts of points of T*dL \ {0}: those who are not in Eq 
and those who are in Eq depending whether R{y, r/) < or not. In the case where 
p e Eqi there exists at most two points of {r^ — |^p} = {77^ = ^(y, ?/')} which are 
in 7r~^({p}), they correspond to the two roots of the equation 77^ — R{yn,y' ,1]'). 
Consider pGT*9L\{0}. 

• If p ^ Eqj one says that p is elliptic. 

• If Card {tt-^{{p}) n {t^ = j^p}) = 2, p is said to be hyperbolic. 

• If Card {t^^^Hp}) n {t^ = |Cp}) = 1, p is said to be glancing. 

We denote by H (rcsp. G) the hyperbolic (resp. glancing) points of dL. We say 
that p e ^ is 

• non strictly glancing if dy^r{p) > 0, 

• strictly glancing if dy^r{p) < 0, 

• diffr active if dy^r{p) > 0, 

• glancing of order k if 

Kiy',o,r,')(9y.^\y.^=o)iP) = 0, < J < fc - 2 and H^[ylo^^,){dy,r^y,=o)ip) + 0. 

We denote by Qk (resp. Qd) the set of points which arc glancing of order fc (resp. 
diffractive) . The assumption that il has no contact of infinite order with its tan- 
gents ((1) in Assumptions [2?6|) consists in assuming 

(A.2) g= U Qk. 

The generalized bicharacteristic are defined by taking the Hamiltonian trajectory 
of Po inside fi and by specifying how the connection is made on the boundary. 
The only problematic points are the glancing ones where the trajectory arrives 
tangentially to the boundary. Indeed, for p G ^, we have y^ — 'i]d — T^iUiV') — 0- 
Recall that the geodesic trajectories are generated by the Hamiltonian fiow Hp^ 
which writes in coordinates (y, rj) 

Hpo = {-^r^'r{y,-n'),2rid,Vy'r{y' ,T^'),dy^r{y,-q')) . 

If P G Qd, dy^r{y, rj') > 0, therefore 77^ and thus yd increase on the trajectory passing 
in p and this trajectory remains in fl. On the contrary, if p ^ Gd, the coordinates 
on drj^ is nonpositive and the trajectory will leave il: if 7]^ decreases, it becomes 
non positive and so for y^- To overcome this difficulty, one uses the vector 

H{p) = Hp„{p) ~ dy^r{y, T]')dna- 

This vector has a coordinate on 9^^ which is 0. One then defines the generalized 
bicharacteristic as follows (see [19] or [3]). 

Definition 1. A generalized bicharacteristic is a continuous map 7 : R ^- Tj*L 
such that there exists a set I of isolated points with 
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• 7(s) e T*L ug for s^I and j{s) eU for sel. 

• for s ^ I , J is differentiable with 

j{s) = Hpjj{s)) if -f{s)eT*Lugd, 

j{s) = H{p) if jis)eg\gd. 

It is proved in [W that these definitions are intrinsic and that if assumption (jA.2p 
is satisfied, then for pa € Tj*LnSo, there exists a unique generahzed bicaracteristic 
such that 7(0) = po- 

A. 2. Propagation near the boundary. Working in space-time variables, one 
first defines tangential symbols by use of the system of local normal geodesic co- 
ordinates: in an open set O where we have such coordinates (j/,??), the function 
'^(2/1 v') G C°°(L X R'') is said to be a symbol of ^™ if a is compactly supported in 
O in the variable y and satisfies: Va, (3 £ N'', BCa^p > 0, 

(A.3) V(y,7y) e I x R^'+i, \dX, {a{y,r^'))\ < C^.p{l + h'l)"-!^! 

This definition implies that the sets A™ depend on the choice of the open sets 
{Oj)j^j. Then, one defines 

A"' ^{aeC°°{TixR'^), Ba'eAT, 3a^eA^, a:^a^ + a'}, 

and one associates with our sequence Un{t,x) its H^ space-time microlocal defect 
measure p by: up to extraction of a subsequence, for all q = q^ + q^ £ A? , 



{q{z,Dz)un,Un) — > / 900(2;, C)dp(^,C) 

where 5'^(2) is the sphere of radius 2 of R'^+i and for (z, () £ 5'^(2), 

qoo{z,C)= lim R^'^q{z,RC). 

fl— S- + 00 

Note that on S*ft, the measure /x is the usual microlocal defect measure (or H- 
measure) as introduced in [12] and [20]. The link between this measure fi and the 
measures /i^ of Theorem 12.41 is described in the following Proposition. 

Proposition A.l. In II x fli x S'^{2), we have 

(A.4) p{t,x,T,0 = -5{t + 1)® p+{x,£,)®dt + -5{t - \)® p-{x,S,)® dt. 

, . 1 {'^{x)uj) 

Hpgp = 4TaaoP- with aoo(x,oj) ^ 



2 B{x)uj ■ uj 

This proposition is proved in section IX3l below. The next step consists in ana- 
lyzing /i on the boundary. We recall that if p G T*dL is an hyperbolic or a glancing 
point, then there exists two points p^ such that 7r(p='=) = p and p^ € {r^ = |^P}. 
These two points are equal if p € 5. When p GH, they differ by their ^-component 
that we will denote by ^^. Besides, as in [4], because of the equation satisfied 
by (un)neN, the sequence of the normal derivatives (i9ArUn)nGN is a uniformly 
bounded family of Lf^^{Il,L'^{dQ)) and we denote by A its microlocal defect mea- 
sure (we suppose that we have extracted a subsequence so that A is uniquely de- 
termined). Following [4], one aims at calculating 

i := Hpgp,- 4TaooM- 
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We already know that if q is an interior symbol, {q, £) — 0. For analyzing the action 
of £ on tangential symbols, one computes the quantity 

In{t) = {[dl+R{y,Dy,),qt{y,Dy,)] Un{t) , u„(i)) 

for qb a tangential symbol of order 1. 

The computation of /„ and the passage n — >■ +00 arises terms on the boundary 
which are the same than in [4] , namely a distribution £1 + £2 where ii is a measure 
absolutely continuous with respect to /i and where, denoting by N(x) the exterior 
normal vector to Oft, £2 satisfies 

''^ (e+-e-)-iv(.x) ''«^^- 

The reader can refer to [3] (pages 14-15) and [4] where the computations are care- 
fully carried out. 

On the other hand, if one uses the equation to transform /„ , there appears a new 
term which was not in the preceding articles. We are going to discuss why these 
terms are harmful when one has (3) of Assumptions 12.61 This term involves the 
quantity V • {'y{x)On). We introduce the principal symbol j{y) ■ t] of the operator F 
which arises when writing V • (7-) in the normal geodesic coordinates and we obtain 

In{t) = {qb{y,Dy,)un{t) , 7(y) • V0„(t)) - {qb{y,Dy,)^{y) ■ V0„(i) , M„(i)) + o(l). 

Let us focus on the first term of /„. In order to study its contribution on the 
boundary, we use S > and a function x G C§°(R) such that x{t) = 1 for i < 1/2 
and x(i) = for t > 1 with < x < 1 and we study 

Jn,s{t) = [x (y) qbiy,Dy,)un , 7(y) • VOnj ■ 

The fact that 7 is tangent to the boundary implies that 

7(y',0)-(7'(y',0),0); 

therefore, we can write 7 = (7', ydld)- Then, the worst term to estimate - which is 
the one which involves dy^ derivatives of 0„ - writes 

JnA'^) == ^ (X (y) <lb{y-.Dy')ydld{y)dy^en(t) , M„(i)j 

and we observe 

< 5\\Un{t)\\Hi(Q:)\\dyen{t)\\L'^i^n)- 

Therefore, for any 6 G C5^(R), there exists a constant C 

limsup / Q{t)Jn,5{t)dt 

n— S-+00 J 



Jn,s{t) 



<C5. 



Finally, letting 5 go to 0, we obtain that this term has no contribution on the 
boundary. As a conclusion, we obtain i! = £1 -|- £2 and we can conclude like in [4] . 
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A. 3. The link between the measure /i and the measures fif. In this section, 
we prove Proposition lA.il 

Note first that the second assertion is a simple consequence of the first one. Suppose 
that we have (jA.4[) where, by Proposition I4.2[ the measures /ij satisfy 

We obtain 

= {-^t^^t - ^ ■ ^xi4) ® S{t + i)(g)dt 

+ {^t^it ~ e ■ ^x^h) ® Sir - 1) (g) dt 

= -2a^fj,+ {x,^)<^S{T + l)(^dt + 2a^fj.^{x,^)(^S{T- 1) «> di 

Let us now prove (|A.4p . We take q G A^ and apply q{z,Dz) to the first equation 
of dH]). We get 

= (g(z, D,) {d^Un - Au„ + V • {l{x)en)) , Un) . 

Using that (0„)„gN is u.b. in Lf^^ (R, iJ^(r2)) and that u„ goes to weakly in 
H^{il), we obtain passing to the limit 

q{z,C){\e~r')dKzX)=0. 
/Rxnix5<*(2) 

Therefore, on the support of ^, we have r^ = |^p. Since moreover 

\C\'^T^ + \e = 2 on 5^(2), 

we obtain that r^ — |^p = 1 on the support of /i, whence the existence of two 
measures fi± on H x ili x S'^^^ such that 

(A.5) 

fi{t,x,T,^) =6{T + l)(E)fi+{t,x,^)+6{T-l)(g)fi^{t,x,^) in V {R x ili x S'^{2)) . 

Let us now link the measure fl± with the measures fi of Section [4.21 inside il. Let 
us consider as in Section l¥^ 

f±,„{t,x) = -j= {idtUn{t,x) ± \Dx\Un{t,x)) . 

These families are uniformly bounded in i^(51) for all i e R and their microlocal 
defect measures are the measures /i . Besides, by the definition of fi, for q G A'^, 
we have 

(?(^,^.)/±,n,/±,„) = ^i{idt±\D\)q{z,D,)iidt±\D\)u„,u„) 
— > If qo.iz,C)i~r±\^\)'d^iiz,0 

«-> + oo 2JmxS''{2) 

= 2/ qooit,x,T^,^)dfl±{t,x,^) 

jRxOxS^'-i 

where we have used (jA.Sp for the last equality. Let us choose q{z,C) = x{'t)K^^O 
with b{x, ^) symbol of order compactly supported in 57, we obtain in the one hand 

/ qoo{t,x,T'i^,^)dfl±{t,x,^) ^ x{t)Kx,Odfl±{t,x,^). 

jRxnxS''-^ jRxnxS''-^ 



ANALYSIS OF THE ENERGY DECAY OF A DEGEN. THERMOELASTICITY SYSTEM 27 

In the other hand, we have for all t G R, 

{b{x,D)f±^n{t),!±,n) — 7> / q{x,uj)(i^if{x,uj). 



Whence fl± (t, x, lo) — -^^ (x, uj) <S) dt. 
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